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Abstract 

Positive pion electroproduction from the deuteron near threshold has been con- 



> 

I sidered within an approach based on the unitary transformation method. The gauge 

O ' independence of the treatment is provided by using an explicitly gauge independent 

expression for the reaction amplitude. The results of calculations for kinematics 

I of the experiments on forward-angle vr"*" meson electroproduction accomplished at 

O I Saclay and Jefferson Laboratory are discussed and compared with those given by 

' the impulse approximation. It is shown that the observed behaviour of the cross 

sections is in accordance with the calculations based on the pion-nucleon dynamics. 

Q I In particular, the pion production rate suppression in the ^II(e, e'7r''')nn reaction 

^ ■ compared to that for the ^H(e,e'7r+)n one can be due to such "nuclear medium" 

• • . effects as nucleon motion and binding along with Pauli blocking in the final nn state. 

> ■ 

X 

^ ; 1 Introduction 

The experimental facilities (such as those at Jefferson Laboratory) put into operation 
for the last decade have opened new horizons in studying the structure of nucleons and 
nuclei. A lot of valuable information can be obtained, in particular, in high precision ex- 
periments on pion electroproduction, including measurements of polarization observables 
and the separation of the structure functions (SF's). These possibilities have already 
triggered a number of novel theoretical approaches (see, e.g., Ref. P) developed for the 
treatment of pion photo(electro-)production off the nucleon. Furthermore, the data [2] on 
neutral pion production on the proton near threshold also renewed the interest in this re- 
action. Applications [3] of the chiral perturbation theory to the threshold pion photo- and 
electroproduction resulted in new conclusions about the predictive force of the low-energy 
theorems, which, in turn, have been rederived 

More understanding of mechanisms of pion photo(electro-)production can be gained 
in studies of the reaction on bound systems of nucleons. Apart from offering means to 
extract the amplitude of pion production off the neutron, the processes like (see Fig. H^) 
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Figure 1: Pion electroproduction on the deuteron: a) kinematics and b) the lA ansatz. 



could shed hght on the details of the nuclear structure. 

For the description of photomeson processes on nuclei, one usually addresses the im- 
pulse approximation (lA) ^ (see Fig. ^p). Within this approach, the two-body mech- 
anisms of the reaction (1) are disregarded along with binding effects in the one-body 
contribution. The latter assumes the validity of the on-energy-shell condition 



for the elementary 7N —>■ vrN' amplitude (the notations are given in Fig. Q). This con- 
dition helps one to provide the gauge independence (Gl) of calculations within the lA. 
However, Eq. (2) may contradict the energy conservation law for the reaction (see, e.g., 
the discussion in Ref. jO]). In general, there are no enough grounds to use the lA in 
studies of the reaction, when, in particular, high-momentum components of nuclear wave 
functions are probed. Moreover, for exclusive and semiexclusive reactions like ^H(7, 7r)NN 
and 2H(e,e'7r)NN with the given pion 3-momentum, Eq. (2) may overdetermine the reac- 
tion kinematics providing a different value for this vector. This leads to violation of the 
transformation properties of observables for these reactions in calculations based on the 



We present in this paper some applications of an alternative formalism developed in 
Refs. El for the description of many-body currents and off-shell effects in the theory of 
photomeson processes on nuclei. Our approach is based upon the unitary transformation 
method (H] in combination with the idea [TIH ITU IT^ to provide the GI of calculations 
through an extension of the Siegert theorem. 

The formalism employed in our consideration is outlined briefly in Sec. 2. The practical 
recipe for providing the GI of calculations is given in Sec. 3. In Sec. 4, our results for the 
cross section of the ^H(e, e'7r"'")nn reaction are presented and discussed. Sec. 5 contains 
concluding remarks. 
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lA (see Ref. 0). 
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2 Underlying formalism 



The present consideration refers to the d(e,e'7r"'")nn reaction. Its amphtude to the one- 
photon exchange approximation (OPEA) is given by 

= [2{2'nfuj]--^ sM) J^M) , (3) 

Jr/g) = (7r+nn;out| J^(0)|d), 

where e^{q) is the (virtual) photon polarization vector [q ■ e{q) = ujeo{q) — qs{q) = 0], 
q = {uj,q) is the 4-momentum transfer (the photon momentum [see Fig. ^]), and J'^(O) 
is the electromagnetic (e.m.) current density at the space-time point x = {t,x) = 0. The 
corresponding GI requirement reads: 

Jtfiq) = . (4) 

Using the reduction technique prescriptions (see, e.g., Ref. one has 

(7r+nn;out | J,(0) | d) = {2{2'KfE^y"^ (nn; out | T,(pJ | d) , (5) 

Up.) = ^ J e^"-^(n + ml) T{<p{x)JM} d'x , 

where 0(x) is the pion field operator. 

Then, following Refs. jHllHI, we apply the unitary transformation method p, to con- 
struct the effective operator acting in the space of states | Xd) (| Xnn)) with no mesons: 

(nn; out | T,(pJ | d) = {xL^ \ [Up^)^ I Xd) ■ (6) 

This operator consists of the one-body contribution T^^l and the contribution T^^l due to 
pion production by a couple of nucleons: 

[Up^)r = Tm + Tl'^ . (7) 

A priori, one cannot exclude the mechanisms such that the photon is absorbed by one 
nucleon, and the pion is emitted by another one. However, the perturbation analysis [S] 
has shown that the corresponding contribution to the lowest degree in the vrNN coupling 
constant g is dropped out, and the perturbation series for T'^l begins with terms of order 
of g^. This cancellation of the simplest two-body pion production contributions is similar 
to the so-called Jennings cancellation mechanism jTHI in vrd scattering. To calculate 
the one-body contribution T^^\ one may use an amplitude of pion production on a free 
nucleon obtained within some dynamical approach making in the corresponding expression 
the off-shell transition — > = (E-^i + Et, — E^, q) (see the discussion in Ref. 0). 
This transition reflects the nuclear medium influence on the one-body pion production 
mechanism. It implies, in particular, that the hadrons participating in pion production on 
a bound nucleon even by a real photon (g^ = 0) are not considered as point-like particles, 
since 7^ in this case. It should be stressed that, unlike the lA (see Eq. (2)), no extra 
kinematical assumptions are made, and the hadron momenta take on values determined 
by kinematics of the reaction on the bound system (see Fig. H^). Thus, the structure 
functions determining, e.g., the differential cross sections of pion electroproduction on the 
deuteron (see Sec. 4) preserve the properties given by their definition (see discussion in 
Ref. [71). 
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3 Gauge independent expression for the amplitude 



In general, calculations based on the contribution only do not satisfy the GI re- 
quirement (4). Further difficulties can stem from the necessity to take into account the 
structure of hadrons, when describing pion electroproduction at ^ 0. Then, an in- 
completeness of the description may lead to results, which are not gauge independent 
even in case of the reaction on a free nucleon (when the lA assumption (2) holds). To 
restore the GI of the treatment, one often adds an extra term to the amplitude making 
the subtraction 

Jf, ^ J^- q- J/q'^ ■ (8) 

Of course, this procedure cannot reflect the complexity of the reaction mechanisms such 
as, e.g., the two-body processes. Moreover, it does not affect the transverse components of 
the transition matrix and is not unambiguous admitting extra subtraction of an arbitrary 
vector such that q ■ X = 0. 

In our consideration, to provide the GI of calculations, we make use of the exten- 
sion [ini HH [12] of the Siegert theorem expressing the amplitude in an explicitly gauge 
independent way through the Fourier transforms of electric {E{q)) and magnetic {H{q)) 
field strengths, 

T,f = E{q)D,, + H{q)M,i , (9) 

E{q) = t[2{27rfuj]--^ {us - eoq) , 
H{q)=z[2{2nfuj]-'- [q x e] , 

with Dif and Mi{ being matrix elements of generalized electric and magnetic dipole 
moments of the hadronic system containing the information on the nuclear dynamics. 
To get representation (9) (see Ref. consider expression 

6{P, + q-Pf){Pf,f\ J^(0) I P„z) = {2nr'l exp(zqs)jf/s)ds , (10) 
j^Js) EE (p,fis),j„{s)) = {PfJ\ J^(s) I P,,^) = {PfJ\ J^(0) I P„z) e-(^/-^-)^ . 



Multiplying the space part of matrix element (10) by an arbitrary vector £{q) and applying 
the Foldy trick yLQj 

Se'l^ = [\vs{ese'^1^) - iXs x [q x e]e'^1^} dX , 
Jo 

with help of the GI condition divjijj(s) = —i{Ef — Ei)pif{s), we get 

6{Pi + q-Pf){PfJ\ J(0) I P„ i) = t{Ef - E,)d,f{q) - t[q X m,f{q)] , (11) 

difiq) = (27r)-3^'dA J e'^'J'sp.fis) ds , 

rriifiq) = (Svr)"'^' AdA J e^^^^ [s x j,/s)] ds . 
Then, due to charge conservation, one may write 

J i{Ef - Ei)qdiM) dP, = u{Pf, f I J°(0) \Pf-q, i) 
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-{Ej - E,{Pf)){Pj, f I J0(0) I Pf, z) = u{Pf, f I J°(0) I Pf - q, t) (12) 

Integration of Eq. (11) over Pi with taking into account relationship (12) resuhs in rep- 
resentation (9), with quantities Dij and Mif being defined as 

In case of reaction (1), one has 

A/ = J\xq[iuj + M^~ ^Ml + (1 - \)W ) J-f{M) ] dA , (13) 

Mij = I f Vxq X Jif{\q) A dA , (14) 

where J^j{\q) is obtained from (7r"^nn; out | 7*^(0) | d ) replacing the deuteron momentum 
by (1 - X)q. 

It should be noted that, whereas quantities dif and rriif in Eq. (11) are singular 
and not proportional to the delta function expressing the momentum conservation, the 
representation (9) is free of singularities. Furthermore, it has been derived here (cf. 
Ref. [TT]) without decomposition of the e.m. current into the part associated with the 
motion of the hadronic system as a whole and the intrinsic current and, therefore, can be 
employed in relativistic calculations. 

This representation generates a correction term additional to the "canonical" expres- 
sion (3), which restores the GI of the amplitude in calculations that fail to satisfy the 
requirement (4). However, when the GI condition (4) does hold, this correction is equal 
to zero automatically. 

In the long- wave limit, Eq. (9) provides the fulfilment of the Siegert theorem ^H] for 
electric transitions in reactions with nonmeson channels ^21- For pion photoproduction 
on the free nucleon at threshold, it leads (as shown in Ref. 0) to the KroU-Ruderman 
result [16^ emerging here as a particular case of the Siegert theorem. The results of the 
application of the approach outlined above to calculation of the differential cross sections 
of forward-angle pion electroproduction on the deuteron are presented in the next section. 



4 Results and discussion 

In the Saclay experiment J7j , determination has been made of the ratio R of the forward- 
angle TT"*" electroproduction cross section for the deuteron to that for the proton at two 
values (1160 MeV and 1230 MeV) of the (virtual) photon-nucleon invariant mass W . The 
incident beam energy E and the electron-scattering angle were 645 MeV and 9 = 36°, 
respectively, while pions were detected at qp^ = 4°. In particular, for W = 1160 MeV 
{u = 290 MeV, | q |=414 MeV/c), the quantity 

d^a-p/dQ'dE'dQ^ ^ ' 

was obtained, with the d^cxd deuteron cross section being integrated over the 50 MeV wide 
peak region. It was claimed in Ref. ^HI that the suppression of pion production rate for 
the deuteron compared to that for the proton is mostly due to the final-state interaction 
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(FSI) between the two outgoing neutrons. However, the corresponding calculations were 
performed for point-like hadrons, while the neglect of the e.m. hadron structure is not 
quite adequate to the kinematical conditions of Ref . [T7] . Particularly, the calculated value 
of 153 pb/MeV-sr^ for the ^H(e, e'7r+)n cross section at W =1160 MeV turns out to be 
much larger than the measured quantity 46 ±3 pb/MeV-sr^. Furthermore, no comparison 
was made in Ref. ^H] between the plane-wave approximation results and those with the 
nn FSI included. 

In the present paper, we apply the formalism outlined above to calculations of the pion 
production cross sections near threshold, in order to estimate the role of such "nuclear 
medium" effects as Fermi motion, binding, FSI of the nn pair, etc. in the quenching of 
the TT"'" meson yield in the ^H(e, e'7r+)nn reaction observed in Ref. |17j . 

The differential cross section of the ^H(e, e'7r+)nn reaction with unpolarized particles 
in the OPEA can be expressed ^19j through the four SF's Wa (« =C, T, 1 and S): 

i4 ^ 



dn' dE' dn^ dE^ . , , V . - ' ^ 2 

+ ^i-^ + v)'^ cos ipWi - |cos2¥;Vrs), (16) 

e = qVq' , V = tan2^/2 , 

where ctm is the Mott cross section, = (cu, q) = {E — E' , k — k'), = cj^ — q^, and 
the SF's are related to the components of the hadronic tensor Wfj^^, by 

WC = Woo , Wt = W,^ + Wyy , Wj = Wo. + W,0 , WS = W,^ - Wyy , (17) 

W,, = {27: fjf 6{q + Pd-p.- Pnn) E;E, (18) 

rm 

with 

Fo = iqDif , 
F = iooDif — iq X Mij . 

One may also write 

F- = J-+G% (19) 
where G" is the correction restoring the Gl, as explained in Sec. 3, defined as 



= ' 4(Ag = 0) , (20) 



UJ 



G 



Vxq{{uJ + Md - v'M2 + (l-A)V)4(Aq) - XqJ^f{Xq)} dX . (21) 

Here we use the reference frame with the OZ axis directed along q, while the momenta 
of incident (k) and scattered {k') electrons define the scattering plane XOZ (see Fig. 
The SF's in Eq. (16) are taken for the angle = between vectors q x and ey. 
Averaging over target polarization states in formula (17) is implied. 

To construct the operator of pion production on a bound nucleon, we have exploited, 
as an example, the so-called generalized Born approximation based on the pseudovector 
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(PV) ttNN coupling (see, e.g., Ref. (201) with the e.m. vertices taken for the on- mass-shell 
hadrons. The corresponding Feynman graphs are shown in Fig. El Parametrizations |2I1 
for the nucleon e.m. form factors were used, while expressions for the nucleon axial and 
pion form factors were taken from Ref. [221 • Unlike calculations ^Hl, the Fermi- motion 
effects have been included in the operator without any nonrelativistic reduction. The FSI 
of the pion as well as the A33 resonance production were neglected in present calculations. 

Using this operator and the explicitly gauge independent expression (9) for the ampli- 
tude, we have calculated the differential cross sections of ^H(e, e'7r+)n and ^H(e, e'7r"'")nn 
reactions for the kinematics of Ref. ^\ at W =1160 MeV. 

For d^crp/df2'd£"df2^, the value of 65.3 pb/MeV-sr^ has been obtained. It turns out 
also that making the gauge subtraction (8) in the amplitude instead of application of 



W=1 .1 75 GeV, = 




0.1 0.2 0.5 0.4 

— (GeV/ c)^ 

Figure 4: Longitudinal and transverse cross sections of the forward-angle pion electro- 
production on the proton in the cm. frame. The dashed and solid curves display results 
obtained with subtraction (8) and representation (9) for the amplitude, respectively. The 
experimental data are from Ref. [221 • 

Eq. (9) leads to the same result, which is some greater than the observed quantity. To 
clarify the origin of this discrepancy, we have calculated the cm. longitudinal and trans- 
verse cross sections of tt~^ meson electroproduction on the proton for the kinematics of 
experiment [2S1 (see Fig. 0}. The results for the longitudinal cross section are in fair 
agreement with data. The description of the transverse cross section as well as our re- 
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suits for the kinematics [T7j could be improved by taking into account the A33 resonance 
contribution to the amplitude, which would affect the transverse component of the pion 
production operator with practically no changes on the longitudinal one. 

When calculating the differential cross sections for the deuteron, we made use of the 
parametrizations [21] for the S- and D-components of the deuteron wave function. The 
nn FSI in the ^So state, which gives a predominant contribution at very small relative 
momenta of the outgoing neutrons, was calculated for the Paris potential exploiting 
version j26| of the matrix inversion method [27j. The results obtained within the approach 
outlined in Sees. 2 and 3 were compared with those given by the lA. In the latter case, the 




Figure 5: Differential cross sections of the ^H(7, 7r+)nn reaction for u = 170 MeV, = 36° 
at low energy of the nn pair. The solid line shows the result obtained within our approach, 
while the dashed (dotted) curve refers to the lA calculations with (without) the ^So-state 
nn FSI included. The histogram exhibits data |28j . 

on-shell extrapolation for the one-body contribution to the amplitude has been chosen in 
a way that provides the fulfilment of the GI requirement (4) at the photon point = 0, 
while, for 7^ 0, subtraction (8) was made. 

In order to estimate the quality of our treatment of the nn FSI, we applied these 
approaches to reproduce Mainz data |2H] on tt"*" meson photoproduction from the deuteron 
at small energies Tnn of the residual nn system (see Fig. It is seen that both the 
approaches give results in reasonable agreement with the experiment. It turns out also 
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that the lA curve practically coincides with the calculations presented in Ref. j28j . 

The plots of the ^H(e, e'7r+)nn cross section obtained for the kinematics of the Saclay 
experiment \Xll are displayed in Fig. El The cross section is presented as a function of 
the "missing" mass (invariant mass of the nn pair). The shown results suggest that, 
contrary to the conclusions of Ref. ^H|; the nn FSI leads to some increase of the 7r+ 



^ 5 



^2 1 
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(MeV) 



Figure 6: Differential cross sections of the ^H(e, e'7r+)nn reaction for the kinematics of 
experiment Notations for the curves are the same, as in Fig. 



production rate within the pion energy range covered in the experiment. Indeed, the 
ratio R (see Eq. (15)) calculated in the plane-wave lA is 0.74, whereas the corresponding 
value given by the lA with the ^5*0 nn FSI included turns out to be 0.79. It is seen also 
that the result obtained using the unitary transformation method in conjunction with the 
extension of the Siegert theorem does not much differ from the one provided by the lA 
indicating a weak sensitivity of the ^H(e, e'7r"'")nn differential cross section to the two-body 
contributions for the kinematics of Ref. J7]. In this case, R = 0.80 was derived being in 
agreement with the measured quantity (15). 

Thus, our calculations suggest that the main source of the 7r+ meson yield quenching 
observed in the experiment ^7] in case of the deuteron target originates from nuclear 
medium effects influencing the 7r+ production process. These effects include nucleon 
Fermi motion and binding along with Pauli blocking of the final nn state. The Fermi- 
motion and binding effects are involved in the one-body pion production operator as well 
as in the phase-space factor arising due to the energy-momentum conservation law, when 
performing integration in Eq. (18). The binding effects in the operator emerge as an off- 
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shell behavior of the one-body amplitude. Their importance in description of exclusive 
pion production off the deuteron was pointed out in Ref. [7|. However, the corresponding 
contribution to the ^H(e, e'7r"'')nn cross section for the kinematics of Ref. J7] appears to 
be not significant. 

The other factors of nuclear medium influence mentioned above prove to be important 
and result in the overall pion production rate reduction by more than 20%. In particular, 
it turns out that the Pauli suppression in the final two-neutron state reduces the ratio 
(15) by about the same value, as the nn FSI increases it. The role of this effect is also 
illustrated by Fig.lZJ where the contributions due to the singlet and triplet final nn states to 
the differential cross section calculated within the plane-wave lA are displayed. It is seen 




Figure 7: Contributions of the singlet (solid line) and triplet (dashed curve) final nn states 
in the ^H(e, e'7r"'")nn differential cross section calculated within the plane-wave lA for the 
kinematics of experiment |17j . 



that the Pauli antisymmetrization of the two-nucleon wave function strongly suppresses 
the triplet contribution in the region of small relative momenta of the outgoing neutrons. 
It should be noted also that the incompleteness of the pion energy range coverage in the 
integral of Eq. (15) also contributes to reduction of the ratio given by that formula. 

For a quasiparallel kinematics and a relatively small momentum transferred to the 
hadronic system (as it takes place in the experiment [IZj), the transverse components T^j^ 
of the pion production operator are dominated by the spin-flip KroU-Ruderman contribu- 
tion jni; whereas the "longitudinal" virtual photons couple mostly to a term proportional 
to the component cr^ of the Pauli spin matrix cr giving no rise to spin-flip transitions. 
Therefore, the singlet nn final state contribution, e.g., in the purely transverse SF's 
and Ws results mainly from the deuteron states with spin projections Md = ±1, while 
the longitudinal hadronic response in the ^5*0 nn FSI peak region is determined by the 
Md = deuteron state. This consideration complies with the conclusions of Ref. JH] 
indicating a strong tt^ production rate dependence on the deuteron spin orientation in 
the vicinity of the ^5*0 peak which, in principle, could be observed experimentally with 
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a polarized deuteron target. As pointed out in paper [TH], this result is predetermined 
by the angular momentum conservation law in the reaction in case of an exactly parallel 
kinematics. We see, that the spin correlations discussed here should take place also for 
moderate deviations from the parallel geometry. However, they might be veiled then by 
the pion FSI and by extra terms in the pion production operator for kinematics probing 
high-momentum components of the deuteron wave function. 

When comparing our numerical results with experiment J7j , it should be kept in mind 
that the data plotted in Fig. El have not been corrected for radiative processes. Introduc- 
tion of such corrections could shift the experimental points by several MeV towards lower 
missing masses making the agreement with our calculations more fair. Radiation broad- 
ening along with finite spectrometer resolution might also be the reason of the fact that 
the FSI peak was not clearly resolved in the experiment. 

Recent measurements of cross sections of forward-angle vr"*" electroproduction on the 
lightest nuclei with the physics motivation similar to that of paper J7j have been carried 
out in Jefferson Lab experiment E91003 j2H|. Our calculations for the ^H(e, e'7r"'")nn 
reaction at one of the E91003 kinematical sets (g^ =-0.4 (GeV/c)^, Vr=1160 MeV, E=8A5 
MeV, and uj= 450 MeV) are presented in Fig. |H1 Near the quasi-free peak, our approach 
leads to the results almost indistinguishable from the ones given by the lA. However, 
we observe an appreciable discrepancy for lower missing masses (greater pion energies) 
that reaches 20 % in the ^5*0 nn FSI peak region. In general, the results of calculations 
are in a satisfactory agreement with the experiment, though in the vicinity of the quasi- 
free peak, the curves are about 15 % above the data. It could be explained by a more 
considerable relative contribution of the transverse SF to the cross section for the 
selected kinematics compared to kinematics of Ref. [T7j. (The ratio (— | + r/)/^^ [see 
Eq. (16)] is equal to 1.7 and 1.4, respectively). Since this SF is sensitive to the A33 
resonance production mechanism neglected in the present calculations, it could result in 
some overestimate of the cross section (cf. Fig. |^. As shown in Ref. [SU], for kinematics 
tuned to the A33 resonance, a considerable cancellation occurs between the resonance 
contributions to the one-body amplitude and the pion rescattering terms. However, for 
the considered above kinematic conditions which are away from the A33 peak for the 
system of the virtual photon and the proton, the A33 contribution to the cross section 
can be still significant through the interference with the nonresonance terms displayed in 
Fig. 121 So, to make a more definite conclusion about the true role of the A33 resonance 
production, a thorough study of FSI including 3-body dynamics is needed (see, e.g., 
Ref. |3T]). We would like to note also that the choice of kinematics for experiments [T7| 0U] 
was predetermined by their physics motivation with a hope to have the cross sections 
strongly dominated by the pion-pole pion production mechanism. It turns out, however, 
that other processes shown in Fig. |21 also play an essential role, in particular, via their 
interference with the pion-pole diagram. In this connection, our calculations confirm the 
conclusions of Refs. [IHl 

5 Conclusion 

The approach based on the unitary transformation method has been applied to the treat- 
ment of pion electroproduction on the deuteron near threshold. To provide the GI of the 
results, we have used the expression for the amplitude stemming from the extension of the 
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Figure 8: Differential cross sections of the ^H(e, e'7r+)nn reaction for the kinematics of 
experiment Notations for the curves are the same, as in Figs. and El 



Siegert theorem. Within the approach, the differential cross section of the ^H(e, e'7r^)nn 
reaction for the kinematics of the Saclay experiment has been calculated and compared 
with that given by the lA. The results suggest that the observed suppression of the vr"*" 
meson production rate in the ^H(e, e'7r"'")nn reaction compared to the case of pion elec- 
troproduction on the proton can be due to such effects as nucleon Fermi motion and 
binding. Calculations performed for the kinematics of Jefferson Lab experiment E91003 
also manifest a satisfactory agreement with the data. Further improvement of the model 
(in particular, its extension to the A33 resonance region) can be achieved by taking into 
account the FSI of the pion with the nn system. 
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